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Some NextGen Control Challenges 




Many redundant effectors 
Surfaces affecting multiple axes 
Actuator rate & position saturation 
Low control authority 
Lighter structures 
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Project Objectives 


Task Objective: Research, develop, and demonstrate a 
flight control system using real-time structural load 
feedback to reduce structural loads on the aircraft and 
remain within its structural limits during normal operation. 



This phase of research will not address gust 
load alleviation, aero-servo elasticity, or flutter 
suppression (March, 2010) 


Subsonic Fixed Wing Project 



Traditional Flight Control System 




Optimal Control Allocation Problem: 

Given B, a desired vector u d and £>0, find u such that 
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Constraints!! 
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Proposed Framework 




Subsonic Fixed Wing Project 
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Optimal Control Allocation 


Optimal Control Allocation Problem: 

Given B, a desired vector u d and c>0, find u such that 
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where F Imax is the load limits at the critical points 


Can also be formulated with loads minimization, with 
constraints given above: 
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load minimization 


Subsonic Fixed Wing Project 
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5.5% GTM-T2 Solid Model 


Subsonic Fixed Wing Project 


Up-scale GTM Simulation 


Up-scale GTM Simulation 

■Simulink model based on 5.5% dynamically 
scaled aircraft derived from wind tunnel & flight 
test data 

■Up-scaled by incorporating Reynolds adjusted 
aero tables 

■Actuator models sized for up-scale GTM 
■NASA Glenn’s CMAPSS (simp2) engine data 
■GTM bare airframe 

■6 ailerons, 4 elevons, 2 rudders, 2 stabs, 2 flaps 
■Vehicle Management System 

- sensor processing module 

- mission manager 

- guidance/control 

■Vehicle Control Augmentation System 

- reference model dynamic inverse controller 

- optimal control allocator 
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Structural Modeling Approach 



Finite Element Model (FEM) 


Assumptions 

■ Static conditions (dynamics will be added later) 

■ Lift is elliptically distributed over wings 

■ Roll is produced by anti-symmetric forces on ailerons 

■ Aileron forces are proportional to surface deflections 

Finite Element Analysis 

■ FEM wing model currently has 90 DOFs (more will be added) 

■ Stiffness matrix K is derived from FEM 

■ K _1 is computed off-line, enabling real-time calculation of 
structural loads 

■ K _1 and B are used to determine T (incremental load matrix) 
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Simulation Test Case 



Altitude 30,000 ft 
Mach 0.85 
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Outline 





Wind Energy Quick Overview 

Utility-Scale Turbines 

Wind Turbine Modeling & Control 

Adaptive Pitch Controller 

Simulation 

Residual Mode Filter 

Results 
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Evolution of Wind Turbines 
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Source: www.owenscorning.com 


■ Wind speed can increase by 20% with a 10 meter increase in height 

■ Average annual wind speed of 6.3 mps (14.1 mph) 50 meters above 
ground needed for economic feasibility 

■ Largest turbine in production is 126 meter diameter (5 MW) 


Wind Power Resources in US 
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Wind Power Classification 


Wind 

Resource 

Wmd Po*er 

Wind Speed* 

Wind Speed* 

Power 

Potential 

Density at 50 m 

at 50 m 
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Fa»* 

300- 400 

64- 7.0 

14.3 - 15.7 



4 

Good 

400- 500 

7 0- 7.5 

15.7 - 16,8 
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Excellent 

500- 600 

7.5- 8.0 

16.8 - 17.9 
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Outstanding 

600 - 800 

8 0- 8 8 

17 9 - 19.7 
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Superb 

800-1600 

8.8- 11.1 

19.7 - 24.8 

“ vv.nd s poods are baaed on a We»ouU k value of 2.0 



• Class 4 or higher wind suitable for 
utility-scale turbines 

• Class 3 areas could have higher 
wind power at 80 meters 




Scary Green Energy Jobs 



Wind Energy Research Areas 
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Small Wind 
Turbines 
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■ Modeling, Design, Materials 

■ Control, Optimization, Acoustics 

■ Sensors, Actuators, Installation 

■ Component Test, Health Monitoring 





Research Areas of Interest to NASA Ames 


Need improved models for wind 
turbine & operating environment 

Wind turbine models do not 
accurately predict performance 
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NREL Wind Turbine at Ames Wind Tunnel 



■ Wind turbines still fail and exact 
cause usually unknown 

■ Mitigation of loads on drive 
train, blades, and tower 

■ Turbine health monitoring 

■ Enhanced energy capture 

■ Design & Optimization 





Utility-Scale Horizontal Axis Wind Turbine 



Utility-Scale HAWT’s 

■ Rotor Diameter: 

- 40-95 m Onshore 

- 90-114 m Offshore 

■ Tower: 25-180 meters 

■ Capacity: 

- 0.1-3 MW Onshore 

- 3-6 MW Offshore 

■ Start up wind speed: 4-5 
mps 

■ Max operating wind speed 
~16 mps 

■ Low speed shaft: 30-60 
RPM 

■ High speed shaft: 
1000-1800 RPM 


Image: NWTC 
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Amount of power 
produced by a wind 
turbine: 

P-^pAC p at 

p = air density 
A = rotor area 
C p = power coefficient 

a) = wind velocity 


Wind Power 


Plot of Power Coefficient versus Tip- 
Speed Ratio 
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= tip - speed ratio 


co 

R = rotor radius 


co T = rotor speed 


Most Basic Governing Equation 


I t cb T =T -T 

rot T aero gen 

T m ,=(\l2)pARC q a,p)or 
TSR : A = Rco, /co 

T ae ro = aerodynamic torque 
< T gen m generator torque 

I rot = rotor inertia 


co T =rotor speed co=windspeed 
p=air density A=rotor swept area 

R=rotor radius X=tip speed ratio 

|3=blade pitch angle 
C q =coefficient of aerodynamic torque 
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Rotation 


source: AWEA 


Wind Flow 


Principles of Wind Turbine Aerodynamic Lift 


Wind Field 
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Aerodynamics 


Dynamic Subsystems 



Internal Loads 


Blade 

Loads 


A 


Structural 

Dynamics 


Motion 

Actuation 



Sensing 


Control 


Electrical 

Power 





Operating Regions & Control Strategies 


Control Objectives: 

Reduce cost of wind energy 
Enhanced power capture 
Mitigate loads to decrease fatigue 
and failure 

Maintain safe turbine operation 
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Pitch-Controlled Wind Turbines 


Region 2: 

- Control generator torque to 
yield optimum power 

- Hold blade pitch constant 

Region 3: 

- Control blade pitch to maintain 
constant rotor speed 

- Generator torque held constant 
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Wind Turbine Actuation 



Blade 

-Pitch 

Nacelle ' 

Yaw* 

/ 

Generator 

Torque 


Control Actions 



Pitch Actuators 



Sensors for Control and Modeling 




• Accelerometer 



Anemometers and 
wind vanes 







Wind Turbine Control and Adaptive Control 


Why is 
control 
important? 


What is 

adaptive 

control? 


Future trends in wind turbines 

- Large multimegawatt turbines placed on tall towers or offshore 

- Softer more flexible systems 

- Increased likelihood of destructive structural dynamic loads and 
responses, excited by highly turbulent flow 

Control can increase efficiency, uptime, and lifespan of turbines 


Uses plant output to modify control law to respond to unmodeled 
plant dynamics, uncertain operating environment and time varying 
parameters 


Benefits of 

adaptive 

control 


Good for poorly modeled plants with uncertain operating 
environments 

Requires less modeling of turbine and its operating conditions 
Reduces controller design time since less tuning is required 
Changes to turbine may require no controller modification 
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Direct Adaptive Control with Rejection of Disturbances 


Plant is linear time-invariant (LTI) 



x = Ax + Bu + ru D 
y = Cx;x(0) = x 0 


where x is plant state, u is control input, y is sensor output, and u D is disturbance input 
and plant parameters (A,B,C,r) are unknown 


Disturbance input comes from a Disturbance Generator: 



z D — IA , z D (0) — z 0 


Disturbances are of known form (i.e.,<fc are basis functions of disturbance), but 
unknown amplitude 

Control objective: Cause Plant output y to asymptotically track 0 

Output error: e, = y-0 — > 0 


Adaptive control Law: u = Ge + G D <p D 


Adaptive gains: G« = -e,e'7? M ,G D = -e y <plh 22 , with h n , h 22 > 0 


For Closed-Loop Stability Analysis, see: Frost, Balas, Wright, IJRNC (2009) 
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Closed-Loop Stability Result 


Theorem : Suppose the following are true: 

1 . All u m are bounded (i.e., all eigenvalues of F m are in the closed 
left-half plane and any eigenvalues on the jw-axis are simple); 

2. The reference model, (A ,B ,C ), is stable; 

5 V m 5 m 5 m 5 

3. Td is bounded (i.e., all eigenvalues of F are in the closed left-half 
plane and any eigenvalues on the jw-axis are simple); 

4. (A,B,C) is Almost Strict Positive Real (ASPR) (i.e. CB > 0 and 
the open-loop transfer function is minimum phase) 

Then the adaptive gains G„, G m , G c , and G D are bounded, and 
asymptotic tracking occurs, i.e. e y = y - y m = Ce* — — > 0 
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Adaptive Pitch Control in Region 3 


■ Objective: Regulate generator speed in Region 3 and reject step disturbances 


■ Controller designed with extended DAC approach to collectively pitch blades and 
hold generator torque constant 


■ Uniform disturbance of wind gust across rotor can be modeled by a step function of 
unknown amplitude, so (t> D = 1 


■ Adaptive control law: 


■ 


u = G e e y +G d 
< = -e y e T y h u 


■ Adaptive controller gains were tuned to minimize generator speed error while 
keeping the blade pitch rate in a range similar to baseline PI controller’s 


Gains used in the adaptive controller were h n = 14.0 and h 22 = 4.0 
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NREL’s NWTC Test Facility 



National Renewable Energy Laboratory’s (NREL) 
National Wind Technology Center (NWTC) Test Facility 

Golden, Colorado 


image credit: 
NREL 


28 


Controls Advanced Research Turbine (CART) 






CART Specifications 


CART Specifications 

■ Variable-speed, two-bladed, teetered, 
upwind, active-yaw 

■ Rotor Diameter: 43.3 m 

■ Hub Height: 36.6 m 

■ Rated electrical power: 600 kW at 42 
RPM in region 3 

■ Region 3 Rated generator speed: 1800 
RPM 

■ Power electronics command constant 
generator torque 

■ Blade pitch rate limit: ±18 deg/sec 

■ Baseline PI Pitch Controller 
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FAST Simulator for CART 




image credit: 
NREL 


■ Aeroelastic simulator capable 
loads of HAWTs 


atigue 

II erodynamics 
^Jtructures 
T urbulence 

predicting extreme and fatigue 


■ AeroDyn subroutine package (Windward Engineering) generates 
aerodynamic forces along turbine blades 

■ Turbine modeled as combination of rigid and flexible bodies 

■ Versatile high fidelity simulation of CART with controller included in 
the loop with switches for DOFs, etc. 
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Simulink Model of Region 3 Controller for CART 


Adaptive collective pitch control for the CART model 



Subtract Operating 
Point Measurements 


FAST Simulation Conditions and Results 



■ Simulation time: 0-100 seconds 

■ Integration step size: 0.006 seconds 

■ Generator DOF switch on, all other DOF switches turned off 

■ Wind turbine had fixed yaw with no yaw control 

■ Aerodynamic forces were calculated during runs 

Step wind inflow Generator speed 
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Generator Speed Error (RPM) 


Results with Additional Flexible Modes Included 



DOF Switches Enabled 

■ generator 

■ drive train rotational 
flexibility 

■ first fore-aft tower 
bending-mode 

■ first side-to-side tower 
bending-mode 


Conclusions 


Adaptive controller was easy and quick to design 
Adaptive controller showed good generator speed regulation 
Pitch rate was comparable for both controllers with step wind 
Adaptive controller performed in robustly with parameter variations 
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Problems with Blade Flapwise Mode 



Baseline 
PI Controller 


Adaptive 

Controller 
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Generator speed with drivetrain torsional and first blade 
flapwise mode DOFs enabled and step wind inflow 


Low-Pass Filter of Generator Speed 



Baseline 
PI Controller 



Generator speed with drivetrain torsional and first flapwise blade DOFs 
enabled with low-pass filter on generator speed 


But low-pass filters have problems ... 
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Residual Mode Filter of Troublesome Modes 


Assume original system: 


into: 


x =A x + B u + F 

p p p P P p D 


V„ =C p x p ;x D ( 0)=x 
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Controlled Plant 


ROM 

All Open-Loop 
Unstable Modes 
that are ASPR 


All Residual Modes 


Q modes 
Open-loop stable 
not ASPR 
in Closed-Loop 


( A,B,C ) is ASPR when CB > 0 and P(s ) = C(sl -A) 'B is mimimum phase. 




Residual Mode Filter Definition and Theorem 


Define the Residual Mode Filter: 



[y Q = 

and the compensated tracking error: 

~ e P = y P -yQ 

Let e Q = *q ~ x q 

Theorem: Let (A,B,C) be ASPR, A Q stable, and 
<p D bounded. Then the Modified Adaptive 
Controller with RMF given above produces 

e v = y n >0 \e n >0 

y y p r— * oo Q t—> °o 

with bounded adaptive gains G e and G D 
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Augmented Adaptive Controller with RMF 


A 



Retained Modes 


x = Ax + Bu 
y = Cx + Du 


NOT ASPR 



Adaptive 

Controller 


Residual Mode 
Filter 



Results with Residual Mode Filter 


Baseline 
PI Controller 


Adaptive 

Controller 




Generator speed with drivetrain torsional and first blade 
flapwise mode DOFs enabled and step wind inflow. 
RMF was designed for drivetrain torsional mode. 
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